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A(X) 
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E 

Q 

p 4 

p” 
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mnxmn matrix 

Xl-A 

lambda matrix or matrix polynomial 
Jordan form matrix 
pseudo-Jordan block 
nxn matrix of Jordan blocks 

primary eigenprojector or matrix residue of [A(X )] _1 
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j>0 secondary latent projector 
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number of generalized eigenvectors for a repeated eigenvalue X^ 
sign matrix 

mnXmn matrix with sign{diag[0 0 ... abs(X^) ... 0]} 
mnxmn matrix of superdiagonal elements for J^ 
ith solvent 

diagonal matrix with ±1 on diagonals 
right eigenvector matrix 

positive projector equal to sum of P^ with Re(X^)>0 

negative projector equ&l to sum of P ^ with Re(X^)<0 

identity matrix 

right eigenvector for X^ 

left eigenvector for X^ 

iii 


T 

P 

adj 
Badj 
det 
Tr 
• ><• 

<• •> 
s 

Q(t,0) 


right latent vector for X^ 

left latent vector for X^ 

transformation matrix 

scalar parameter 

adjoint of matrix 

block adjoint of block matrix 

determinant of matrix 

trace of matrix 

outer product 

inner product 

Laplace variable 

state transition matrix 


Abstract 


The mathematical theory for decoupling mth- order matrix differential 
equations is presented. It is shown that the decoupling procedure can be 
developed fro^ the algebraic theory of matrix polynomials. The report 
discusses the role of eigenprojectors and latent projectors in the decoupling 
process and develops the mathematical relationships between eigenvalues, 
eigenvectors, latent roots and latent vectors. It is shown that the eigen- 
vectors of the companion form of a matrix contains the latent vectors as a 
subset. The spectral decomposition of a matrix and the application to 
differential equations is given. 


1. Introduction 


The purpose of the material in this repovt is to formulate the algebraic 
theory of systems and application to spectral decomposition and decoupling 
of differential equations. The relationship between eigenvalues, eigenvectors, 
latent roots and latent vector*, of matrix polynomials will be given. Since 
most of the equations of motion of vibrating systems are cast in second-order 
form, the algebraic properties of second-order matrix polynomials have an 
important role in the determination of solutions of vibration problems. 

Although the mathematical development will be In general form, the analysis 
Includes second-order matrix polynomials, 

The concept of scalar residues is well understood in complex variable 
theory and the Inversion of Laplace transforms. The theory and use of matrix 
residues is not widely used and the relationship to eigenvectors and latent 
vectors has received little attention. Matrix residues, eigenprojectors 
and latent projectors are useful in analyzing matrix polynomials and time 
domain solutions to differential equations . Several papers have been pub- 
lished in recent years on matrix polynomials, see Dennis, Traub and Weber, 

[1], as well as a short paper by Denman, [2]. Some material on matrix resi- 
dues has been given by Zadeh and Desoer, [3], and on projectors, Cullen [4]. 
Lancaster's book, [5], is an excellent source on latent roots and latent 
vectors of matrix polynomials which he denotes as lambda matrices. 

It will be shown that eigenvalues and eigenvectors of tlve matrix com- 
panion form and latent roots and latent vectors of a matrix polynomial are 
related. The matrix residues of the inverse of a matrix polynomial 

MX) = lX m +A, A"*” 1 *. . .+A , which will be called latent projectors, are sub- 
1 m 

matrices of the matrix residues of the inverse of (Al-A ) where A is in 


1 


2 


block companion form. The latter residues will be referred to as eigen-* 
projectors. It will be shown that latent projectors and eigenprojectors are 
useful for solving simultaneous differential equations. 

The theory of Laplace transforms will be useful in introducing the 
concepts that are to follow. The material on Laplace transforms in most 
textbooks is limited to scalar problems and functions which is unfortunate 
since modern engineering problems are likely to be formulated as matrix 
problems due to the complexities of the systems to be analyzed. The exten- 
sion of Laplace transforms to matrix functions is a simple task provided 
that the development of Laplace theory Is based on algebraic functions 
rather than scalar functions. 

Let f(t) be a scalar function for which the one-sided Laplace transform 
is given by 


( 1 . 1 ) 


L[f(t)J - F(s) = 


r 

■ 0 


f(t)e~ st dt 


with the usual assumption that the Integral exists. The inverse transform 
is defined as 


( 1 . 2 ) 


L" l [F(s)] = f(t) = 2~~j 


fc+i°° 


c-i“ 


F(s)e St ds 


where c is properly defined to enclose all singularities of the integrand. 
If F(s) has the property that 

(1.3) ^|F(s)|=0 

the inverse transform of F(s) is 


3 


<1.4) l” l lF(s)] - f(t) " l residues of [F(s)e 8t ] 9 T p._e 1 

i-1 b- 8i 1-1 

where F(s) is a ratio of two scalar polynomials. 

The Laplace transform method is valid for vector and matrix functions 
provided that certain restrictions are satisfied. Let A(s) be the matrix 
polynomial 


(1.5) 


A(s) - Ia m + A 1 s m ” + ... + A 
1 m 


where all coefficients A.^ are n*n. The inverse of A(s) is in general form 


( 1 . 6 ) 


[A(s)]" 1 


adjtA(s)] B(s) 
det[A(s) ] d(s) 


with adj the adjoint and det the determinant of A(s) respectively. The 

characteristic equation of A(s) is given by d(s) * det A(s) and will have 

a maximum of mn roots, Lancaster [5] calls these latent roots. The Inverse 

-1 

transform of F(s) * [A(s)] when the roots are distinct are 


(1.7) 


- i , mn 

L X [F<s) 3 “ f(t) » L i(A(s)] A 3 - l residues of [F(s) 


Sti 


i=l 


s-s, 


which can be expressed as the matrix analog to (1.4) with 

mn A 

(1.8) f(t) » l F 10 ex P( s ^ fc ) 

A 

where the matrices are matrix residues or latent projectors. It is 
obvious that the latent projectors are coefficients of the partial fraction 
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expansion 


(1.9) 


[AC)!' 1 - I 

i-l ‘ l 


The usefulness of the above approach to Laplace transforms can be 
illustrated by considering n simultaneous differential equations of m-th 
order, i.e. 


( 1 . 10 ) 


.m ,m-l 

A_ — ~ + A. — — j + , . . 4* A x « 0 
0 dt n 1 dt"*" 1 n 


x(0) - c 


where x(t) is a nth-order vector. It follows that 


( 1 . 11 ) 


A(s)X(s) “Ac 
m 


x(0) - x(t Q ) «... x w (0) - 0 


( 1 . 12 ) 


X(s) - [A(s) ] Ac, 
m 


The time domain solution to (1.10) is then given by 


(1.13) 


x(t) ‘ >i 5l ° A, " c exp<8lt> 


when the latent roots are distinct. 

The matrix polynomial given in (1.5) can also arise from the canonical 
form or companion matrix. It is not difficult to show that the matrix A(^) 


given by 
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r xi 

-1 

0 

• * « d 

0 


0 

XI 

-I 

• • • 0 

0 


p- 

p 

9 

• • 

• 

(1.14) 

A(X) - 

* 

• 

.• •» *■ 9 

• 


* 

• 

9 

• 99 .#■ 

• 


. 0 

0 

0 

• 9 9 XI 

-I 


1 A 

A m-1 

A «-2 

k 

9 9 • 

Xl+A. 

1 . 


has the same characteristic equation as (1.5) when A Q - X. The root of the 
characteristic equation obtained from det[A(X)] » 0 will be the eigenvalues 
of A and are equal to the latent roots of A(X) when Aq - I. The eigenvectors 
of A must be related to the latent vectors of A(X); that relationship will 
be given later. 

If z(t) Is defined as the vector 

(1.15) [z(t)]' £ - {[x(t)] T [x(t)) T ... tx <ra) (t)] T } 
then z(t) satisfies the equation 

(1.16) z(t) - A z(t) 


with 

(1.17) Ix (m) (t) + A 1 x <m_1 \t) + ... + A m x(t) - 0, 

The solution vector z(t) is given by 

mn 

z < c ) " I p iQ C exp(X i t) 


( 1 . 18 ) 


6 

where P iQ are the matrix residues of (A(X)J which will be called eigen- 
projectors. It Is assumed that ( 1 . 18 ) is for distinct eigenvalues of A. 
Since the eigenvalues of A and the latent roots of A(X) are the same and 
the eigenvectors of A and the latent vectors of A(X) are related, the 

A 

elgenprojectors and the latent projectors P^q must be related. The 
vector C In (1.18) Is obtained from (1.12) and the definition of the 
canonical form for the Bystem. 


2, Eigenprojectors of Matrices 


Lot A be defined as a nr»xon matrix with eigenvalues right eigen- 
vectors y A and left eigenvectors z^. Define Q as a unXmn matrix constructed 
from the eigenvectors such that 


( 2 . 1 ) 


Q- ly x y 2 y 3 ••• y m ) 


where it is assumed that the columns of Q are linearly Independent and spans 
the c mnXmn space. The matrix Q has the property that a similarity trans- 
formation on A with Q will reduce A to the Jordan form 


( 2 . 2 ) 


Q AQ 


with J ■ diag[J^,J2» • . Jp3 with J^ a Jordan block. The Jordan form will be 
diagonal if A has mn distinct eigenvalues or if A has mn linearly independent 
eigenvectors satisfying [A,^I-A]y^ “ 0. If A has repeated eigenvalues and 
is defective, the Jordan blocks leading to the defectiveness of A will have 
one or more plus ones on the superdiagonal on a Jordan block. It will be 
necessary to utilize the chain rule for generating the generalized eigen- 
vectors for the defective Jordan block. 

The mathematical analysis will be simplified if all eigenvalues with 
the same values are considered as a pseudo-Jordan block with the? plus ones 
on the superdiagonal. Assume that J^ and J i+;L are as shown with J i+ ^ having 
the same eigenvalues as J^ but where J^ + ^ has the plus ones on the super- 
diagonal. Although the definition 
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8 


(2.3) 


i+1 


0 

0 

0 


0 

h 

o 

0 


0 

1 

X, 


1 J 


Is not conventional. The two blocks J, and J. will be considered as a 

1 1+1 

pseudo- Jordan block. 

Let A have q values c* X^ with q pseudo-Jordan blocks as defined. Assume 
that mn-k of the eigenvalues are distinct and q-nwvHc are repeated. Each re- 
peated eigenvalues will have multiplicity r^ and the number of generalized 
eigenvectors for the repeated eigenvalues will be Z^, the number of plus ones 
on the superdiagonal of J.. , It will be assumed that the ones are located in 
the last Z ^ rows of the repeated eigenvalue pseudo-Jordan block. The term 
pseudo-Jordan block will be dropped in the following discussion and the term 
Jordan block will be utilized with J denoting a pseudo-Jordan block. 

In addition to the above assumptions, let denote a mnXmn matrix. 

The first subscript denotes the eigenvalue to which the matrix belongs 
and the second subscript is an index which has a maximum value equal to the 
number of generalized eigenvectors required for the eigenvalue X^; this will 
be Z If Z ^ is the number of generalized eigenvectors for X^, then 
j = 0,1,2, ... ,&^ with j = 0 for zero generalized eigenvectors. If X^ is a 
distinct eigenvalue, then F^q will be defined as 


(2.4) 


F^y = diag[0,0, .. .0,1,0. . .0] 


with the one located in the same row and column as X^ is in J, The matrix 
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/ 


F 10 for a re P eated eigenvalue X ± witli multiplicity ^ - 3 will have cite form 

(2.5) F 1() - diagtO, 0,... 0,1,1, 1,0... 0) 

with the one located in the same rows and columns as X^ is in J. 

'l'o complete the definition of assume that X i has multiplicity 
with r^-A^ linearly Independent eigenvectors and generalized eigenvectors • 
The associated Jordan block will have ones on the superdiagonal with the 

ones located in the last rows of The matrix F^q will be as given in 

(2.5) but the set of matrices ^3 ,F 12 > * * ,F l^i now ®xist with having 
only the ones of the superdiagonal of located on the superdiagonal of F^, 
The next matrix in the sequence F^ will be generated by moving the ones on 
the superdiagonal of F^ up one diagonal position by moving to the next 
columns of F^. To illustrate the, construction of , let J be defined as 
in (2.3), then F. f) F ( ( and F^ are 




1 

0 

0 

r 

0 


0 

0 

0 

0 

* 

(2.6) 

f IS 

lO 

0 

0 

1 

0 

0 

1 

0 

0 

F ** 

il 

0 

0 

0 

0 

1 

0 

0 

1 

F m 
12 



0 

0 

0 

1. 


0 

0 

0 

0 . 



0 0 0 0 
0 0 0 1 
0 0 0 0 
0 0 0 0 


The eigenprojectors of A, [6], can now be defined using the established 
notation. Assume that A has mn-k distinct eigenvalues and q-mn+k repeated 
eigenvalues. The primary eigenprojectors will be defined as 


(2.7) J ’ i0 “ (} K f0 Q ~ 1 ■* 1 , 2 , . . . q 

with the primary eigenprojectors having the properties 


( 2 . 8 a) 

3 


) p * r 

, 10 
i -1 


( 2 . 8 b) 

P P m p 

10 10 l 10 


( 2 . 8 c) 

P p * 0 
.10 j 0 u 

i + J 


<! 

These properties follow directly from (2.6) since ][ F (0 « 1, *• 

ami F JA F., A * 0, The primary eigenprojectors are idempoteut matrices, i.e. 


10 JO 




I’iO “ P,^q where a is a positive integer. 

The spectral decomposition for A which is not defective- is given by 


( 2 . 9 ) 


q 

A * I A.P 

i*l 


which follows directly from the definition of the primary eigenprojectors. 
If A is defective, a set of secondary eigenprojectors will be defined 
as the projectors constructed from the eigenvectors and the sequence of f 

***J 

matrices. Let the secondary eigenprojectors he defined by 
(2.1,0) V L] * Q P|j tf 1 j - 1,2,...^ 


wlLh i set by the associated eigenvalue which is repeated but defective. It 
follows from the definition of that 

(2.11a) P i 0 P ij " l ‘ij J - 1,2,...^ 


p p 

ij ij 


i, j+T 


(2.1.1.b) 


j / 0 
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<2. 11c) 


P lj P i,j+l " 0 


j 4 o 


The secondary eigenprojectors are required for the spectral decomposition 
of the most general matrix A. 

If A has general form and i3 defective then 


( 2 . 12 ) 


A s Q J Q -1 - Q[A + f F. ]Q -1 


i=l 


ir J 


but this is equal to 


( 2 . 13 ) 


A - Q( j, 


or finally 


( 2 . 14 ) 


A “ ^ ^loV*!! 1 


1=1 


Although the secondary eigenprojectors with j>l are not necessary for the 

spectral decomposition, it will be shown later that the partial traction ex- 
~1 "I 

pansion of [A (A) ] = [ AX— A ] can be expressed in terms of the eigenprojectors, 

The procedure given for computing the eigenprojectors for a matrix has 
been based on the assumption that the right eigenvector matrix, Q, is known 
completely. The inverse of Q will have row vectors that are the left eigen- 
vector of A thus P^q depends on the right and left eigenvectors, y^ and 
respectively. The right and left eigenvectors for distinct eigenvalues are 
determined from the equations 


( 2 . 15 ) 


t A i X-A]y ± = 0 


\ 


12 


(2.; 16) 


z i [X i l-A] - 0 


or equivalently ** 0 for the left eigenvectors. If A is defective 

for an eigenvalue X^, the chain rules, 17] 


(2.17) 


[X,I-A]yk +1 « -y^ 


k “ 1,2, . . .i. 


(2.18) 


r , T ,T W IcflvT , k N T 

[X^-A ](z i ) - -(z t ) 


k « 1,2, . . 


1 1 

are used for the generalized eigenvectors where y^ and z^ are any one of the 
linearly Independent eigenvectors for the repeated eigenvalues. 

The primary eigenprojectors were defined earlier and are given by 


(2.19) 


P i0 ■ Q F i0 Q 


Let Q = Q be the matrix of right eigenvectors and be the matrix of left 
eigenvectors for distinct eigenvalues with 


( 2 . 20 ) 


Q r - [y L y 2 ••• y^l 


and let be scaled such that Q^Q r = I. It then follows for a distinct 
eigenvalue X. that 


( 2 . 21 ) 


P i0 - Q r P i0 - Vi • 


I- 
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Since Q r F i0 Qjj will be given by the outer product of y^ and ‘z ±f the scaling 
of is equivalent to scaling y^ or z^ such that <2^7^ » z^'y * 1. Since 
any eigenvector can be multiplied by a constant then any set of arbitrary 
scaled eigenvectors can be used to compute provided that the arbitrary 
constants are removed by dividing by the scaling factor. If y^ and 2 ^ are 
arbitrary eigenvectors the elgenprojectors for the distinct eigenvalues are 
given by 


( 2 . 22 ) 


10 


y i 2 l 

Vi 


The elgenprojectors for the repeated eigenvalues when A Is not defective 
are determined from a simple extension of (2.21). Since the eigenvectors 
are linearly independent, is given by 


(2.23) 


10 


Wit 


l y“z 
J-l 


JrJ 

i 2 i 


J j 
r i 

i^i 


l 

d-i 


where the superscript denotes the jth eigenvector belonging to the pseudo- 
Jordan block. 

The elgenprojectors for the repeated eigenvalues when A is defective 
are computed in a similar manner to the repeated eigenvalues for the nonde- 


fective case. The. primary eigenp rejector for a defective Jordan block is 
given by (2.23). The secondary elgenprojectors can be computed from the 

/s 

eigenvectors by considering (2.10). Let F. . denote the subblock of F for 

^•J **-J 

a Jordan block J^ where F^. is r^xr^. Let ^ and ^ denote the rectangular 
matrices of right and left eigenvectors respectively of A^. Equation (2.10) 


can be rewritten as 


14 


r n 


(2.24) 


lj 


Y lVi 


[y'¥ 

l/ i y i 


_ r i * 

^ lF « 


-2 



where Q f and have been properly scaled. Let f gt denote the elements of 

A A 

with f gt « 1 or 0 depending on the 1th row and jth column of F^. Equa- 
tion (2.24) can then be written as 


(2.25) 



¥ 



ij 



I 1 I 1 •• 8t 0 21 ] 

s-1 t=l z® y® 


A 

Since F^ will be sparse, only a few terms of the summation are required. 

An example will now be given to illustrate the computational procedure 
using the eigenprojecters. Let A be given by 



with Jordan form 


10 0 0 
0 3 10 

0 0 3 1 

0 0 0 3 



The right eigenvector matrix Q is 


15 


Q 


r 


1 

1 

1 

1 



r 12 3, 

[y l *2 >2 > , 2 1 


and the left eigenvector matrix by 



The eigenprojector is found from y^ and z^ and Is 


10 


y l Z 1 
2 1 y l 


1 

4 


1111 

1111 

1111 

1111 


10 


wltlt the primary eigenprojector for A2 = 3 given by 


20 


- I 


v J z J 
y 2 '2 


j-l ZJ y J 2 


1 1 
y 2 Z 2 


T + 


Z 2 y 2 


2 2 
y 2 Z 2 
2 2 
Z 2 y 2 


3 3 

y 2 Z 2 
3 3 

*2 y 2 


3 -1 
-1 3 

-1 -1 
-1 -1 


-1 -1 
-1 -1 
3 -1 
-1 3 


F = I 
20 3X3 


Ah a (‘heck, on the two eigcnproj ectors , note that P^+Pgo “ !• The first 
.secondary eigenprojector is given by 


y? f, 


1 2 
y 9 z 9 


2 3 
y-j z o 


y y 7 2 *11 2 '2 "2 , y 2 "2 
A A s t 1 ,1 2 2 
s 1 t 1 z 2 y 2 z 2 y 2 z 2 y 2 


0 0 
2 -2 


2 -2 

0 0 


f 23 " 1 


The second secondary elgenprojector Is determined from (2.25) with f , „ « 1 

X J 

/s 

as the only nonzero element in P 09 . The elgenprojector P 99 is then 


y 2 Z 2 
1 1 
Z 2 y 2 


1 -1 
1 -1 


1 -1 

1 -1 

-1 1 

-1 1 


A check will show that P^g P^^ * P 21 ,P 20 P 22 " P 22 an< * P 21 P 22 “ °* 

The spectral decomposition of A is given by X 1 P 10 + X 2 P 20 + P 21 which 


1111 

1111 

1111 

1111 


3 -1 
■1 3 


10 

-2 

0 

-4 

0 

8 

-2 

-2 

-4 

0 

10 

-2 

-2 

-2 

-4 

12 


which agrees with A as given. 
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The elgenproj actors can also be computed from the inverse of A(A) where 
A(A) - X'C-A. The most general partial fraction expansion of LAMjf 1 is 
given by 


(2.26) 


1 q P .«A K 

i-l _ V / 10 . ri i 


tA(X)f - l <*=r + l Ti+i ] 

1-1 A A i j-1 (A-A.) j+1 


The three cases, distinct eigenvalue, repeated eigenvalues with multiplicity 
r^ but not defective and the defective matrix must be discussed. 

The three cases can be analyzed by considering the Jordan blocks for the 
three different eigenvalue cases. Rather than consider the mixed Jordan 
forms, consider the three individual Jordan forms. Assume first that 
J « diag[X^ X 2 X^] thuB 


(2.27) 


x x ° ° 

A * QJ(f 1 * Q 0 X 2 0 (f 1 

0 0 X„ 


with [A (A)] A given by 


(2.28) 


tA(X )]" 1 - [Xl-Af 1 - Q 




(X-X 2 ) 


(X-X 3 )' 


It follows directly from (A-A^) [A(A)]’*' 1 evaluated at A^ that 


(2.29) 


U-yiU-Af 1 « Q F i0 Q - 1 - F i0 
A-X^ 


Consider the Jordan form with X^ repeated 3 times but where A is not 
defective. The inverse of A (A) is 


r 
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(2.30) 


[AU))" 1 - [Al-Af 1 - Q 


(A-XjT 1 0 0 

0 (A-Aj)" 1 0 

0 0 ' <A-A 1 >” 1 



Since A^ has multiplicity 3 then (A-'A 1 ) 3 [A(A) ] -1 will be 


(2.31) 


(A-A 1 ) 3 [A(A>r 1 - Q 


(A-A.) 2 0 0 

X rt 

0 (A-A.) 0 

0 0 (A~Aj) 2 



It is obvious that (2.31) will be zero when evaluated at A - A^. The first 
derivative of (2.31) with respect to A will also be zero at A - A^ with the 
second derivative of (2,31) given by 


(2.32) 


^ {(x-v 3 [A(X)f l ) 

dA 


Q 


'2 0 
0 2 
0 0 


2Q F l 0 Q 


-1 


The eigenprojector for the repeated eigenvalue matrix with multiplicity r^ 
but for A not defective is given by 


(2.33) 


= Q^qQ" 1 - 


10 


V 1 r 

dijiV tx -v liAwrl) 

1 dA 1 


A-A, 


The defective matrix will be analyzed by considering the Jordan block 
with A^ of multiplicity 3 and 2 generalized eigenvectors. Let [A(A))"* 1 be 
given by 


A-Aj 

0 

0 


-1 

A-A x 

0 



-1 


tf 1 "*} 


A-A x 


(A-A^" 1 

0 

0 


(A-AjT 2 (A-A^ 
(A-A 1 ) -1 (A-A x ) 
0 (A-A 1 ) 


(2.34) 


[A(A)]” 1 - Q 
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The sequence of evaluations of (X-X,j) IA(X)J and the derivatives with re- 
spect at X «* A^ will be 


(2.35) 


(X-X 1 ) 3 [A(X))“ 1 


x*x. 


0 0 1 
0 0 0 
0 0 0 


<f 1 - PP^Q - 1 - P 12 


(2.36) ^ {(X-X 1 ) 3 [A(X)]" 1 } 


x»x. 


(2.37) \ UX-XJ^AU))- 1 } 

1 dk l 1 


A"A, 


0 10 
0 0 1 
0 0 0 

10 0 
0 10 
0 0 1 


<1 ' 1 ■ v 1 ■ p u 


Q' 1 - QP 10 Q' 1 - P 10 


The computational procedure for finding the elgenprojectors, usually 
referred to as matrix residues In the above, from the partial fraction ex- 
pansion can now be summarized. If A.^ is distinct then 


(2.38) 


10 


(A-A^ [A (A)] 


-li 


X«A, 


with the primary elgenprojectors for the repeated eigenvalues A^ of multiplicity 
r^ given by 


(2.39) 


^i"" 3- r 

p i0 - V-t [(X -V 

1 dA i 


A**A, 


The secondary elgenprojectors are defined only for repeated eigenvalues with 
multiplicity r^ with A defective and requiring generalized eigenvectors. 
The secondary elgenprojectors are given by 


(2,40) 




1 d 

J' *1 


j t(X-Xj ri [A(X))' 1 ) 


J * Of 1| i i i ^ 


"1 

where (A(X)J has all common factors of adj[A(X)] and det[A(X)) cancelled * 

—1 A *“4 

so that [A(X)1 Is a minimum polynomial with r^ equal to the power of (X-X^) 
in the denominator. 

The computation of the elgenprojectors by the residue method will be 

— i 

Illustrated with the previous example. The Inverse of [A(X)] Is given by 


£A<X> r 



1 0 o o" 


-7.3 -0.5 0 -1 

1 

0 10 0 

x 3 + 

0 -8 -0.5 -0.5 

(X-l) (X-3) 3 

0 0 10 

-1 0 -7.5 -0.5 


0 0 0 1 


-0.5 -0.5 -1 -7 , 


17.75 

3.25 

1.25 

4.75 


-13.25 - 4.75 

- 3.25 

- 5.75' 

0.75 

20.25 

3.25 

2.75 

X+ 

- 2.75 -15.25 

- 4.75 

- 4.25 

5.25 

0.75 

17.75 

3.25 

- 6.25 - 2.75 

-13.25 

- 4.75 

3.25 

2.75 

4.75 

16.25. 


- 4.75 - 4.25 

- 5.75 

-12.25 _ 


The elgenprojector for X * 1 Is 


10 


1 

4 


111 

111 

111 

111 


(X-l)lA(X)] 


-1 


X»1 


The eigenvalue X * 3 will have a primary elgenprojector and two secondary 
elgenprojectors. The elgenprojectors ? 22 , P 2 ^ and P 2Q are 


21 


P 22 " U-3) J [A(X))* 


P 21 " dX ta -3) 3 [ A < X > 1" X } 


1 d 


P 20 2 2 {(W) [AW] } 


dX 




-1 

1 

1 

-1 


1 

-1 

1 

1 

-1 

X -3 

4 

1 

-1 

“1 

1 



1 

-1 

-1 

1 



* 

0 

0 

2 

-2 


. . 1 

2 

-2 

0 

0 

X -3 

4 

-2 

2 

0 

0 



_ 0 

0 

-2 

2 



m 

3 

-1 

-1 

-1 


« A 

-1 

3 

-1 

-1 

X -3 

4 

-1 

-1 

3 

-1 



-1 

-1 

-1 

3 


The eigenprojectors agree with the previous values round from the eigenvectors, 
The partial fraction expansion of [A(X)) Is 


(2.41) 


[ A ( X )]- 1 


* 




+ 


p p 

21 + 22 

( X ~3) 2 ( X -3) 3 


Two methods of computing the eigenprojectors (or matrix residues) have 

been discussed in this section. The first procedure given was based on the 

•~1 

eigenvectors with the second method requiring the inversion of £ A (X) 3 and 

evaluation of the residues. It has been shown that the two methods are 
equivalent although the numerical computations may not necessarily be com- 
parable. 


3. Lambda Matrices and Latent Projectors 


- 


The analysis in Section 2 was based on the assumption that the A matrix 
was in general form with mn rows and columns. This assumption is valid when 
the differential equation describing the dynamics of a system are in the 
state variable form, first-order differential equations. The system equation 
may not alw* ,f s be in first-order form as it is common practice in some engi- 
neering disciplines to write the differential equations in mth-order form. 

If such is the practice, then lambda matrices will be encountered, This 
section considers lambda matrices or matrix polynomials, Gantmacher, [8]. 

Assume that A(X) is a matrix polynomial in X of mth-order with nxn co- 
efficients of the form 

(3.1) A(X) - A () X m +A 1 X m ' 1 +...+ A^X+A 

which Lancaster, [1], calls a lambda matrix. Dennis, Traub and Weber [3], 
make a distinction between A(X) and A(X) where X is n*n by calling the 
latter a matrix polynomial. The polynomial in (3.1) is commonly referred 
to as a matrix polynomial in control theory and that designation will be 
followed here. 

The roots of det[A(X)] are called latent roots and the vectors that 

— kS 

satisfy [A(X^)]y i are referred to as latent vector. This terminology will 
be followed in this work to avoid confusion with eigenvalues and eigenvectors. 
The concept of latent projectors will be introduced in this section with 
the latent projectors having an analogous role to eigenp rejectors . It will 
be assumed that a latent root may have multiplicity r^ and that A(X) will 
be defective requiring generalized latent vectors where a defective 
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lambda matrix has the same meaning as A being defective, there will not be 
r^ linearly latent vectors for the latent root 

If A q Is Invertible, then (3.1) can be written as 


(3.2) 


A(X) - A 0 [lX m +A l X m " X +...+A n| _ 1 X+A iu ) - A q A(X) 


where A^ *• Aq A^. The discussion that follows will focus on A(X) although 
a complete treatment of lambda matrices should Include the case when Aq Is 
singular. 

'l’he latent roots of A(X) will be denoted by X^ with the right and left 

A A 


latent vectors, 


denoted by y^ 


and z^ respectively. 


The latent vectors for 


the latent roots X^ satisfy 


(3.3a) A ^l^ y i " 0 

(3.3b) Si A(X ;L ) - 0 


for the right and left latent vectors respectively when X^ is distinct or 
A(X^) is not defective. If A(X) is defective for Xi then a chain rule must 
be employed. Lancaster and Webber, [9), have given the chain rule as 

1 2 . * . v i 1 

dX 


(3.4) 


A(x i )y i + ’’dx— y i * 2 - y i + - • * ' 


4j-l 


d~ *A(X.) 
dX 


with y^ as a linearly independent latent vector, 
pendent latent vector than y is given by 


If y^ is a linear inde- 


dA(X^ A 1 


(3. 5a) 


and by 
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(3.5b) 


A < X i> y i 


“dT“ y i 


x d 2 A(X 1 ) 


"1 


with all others computed by recursive use of (3.4). The chain rule for the 

A*| 

left latent vector is similar except that is a premultiplier of the terms 
in (3.4). 

The computation of the latent vectors is described in the example. Let 
A(A) be 


A(A) ■ 

which has latent roots 
Let A * 1 then 


1 

0 

0 

-4.5 

1.5 



5.5 

-3.5 



A“+ 







0 

1 


1.5 

-4.5 



-3.5 

5.5 ^ 




" 

ml 





*1 

“ 1 

X 2 “ 

2, A 3 

m 3 and 

X 4 “ 

X 3 " 

3 with 


1 . 


AU)^ 



with y^ as given for the latent vector. The latent vector for A *2 is ob- 
tained from 


A 

0.5 -0.5 ~ 


0 

A 

-1 

AWy, - 

r 

1 

o 

e- 

o 

t 

y 2 “ 

0 

y 2 " 

-1 


The linear independent latent vector for A * 3 is found from 


A(3)y* 


"l 

1 " 

*1 

o' 

*1 

1 

3, 

1 

y 3 " 

0 

V =3 

y 3 

-1 


and the generalized latent vector for A = 3, computed from the 
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chain rule. Using (3.5a) 


A(3)?3 


dA(3) -1 
dX y 3 


1 1 
1 1 


'3 

'2 


0 

0 


or (y^) 1 - [-1 1]. The latent vectors y 2 and y* are chosen with signs oppo- 
site to y^ and y^ respectively for convenience. Since y^ is in a two-dimen- 
sional space, only two of the latent vectors are necessary to span the space. 
It should be noted that the term linearly Independent latent vector is not 
proper terminology and will be dropped in favor of latent vector hereafter 

a A 

if y^ satisfies MX^)y i ” 0. 

The concept of distinct and repeated latent roots as well as a defective 
matrix polynomial will be clarified by relating A(X) to the companion form 


(3.6) 


A « 


A. It 

- 

is well 

known that 

the 

0 

I 

0 

0 

0 

0 

I 

0 

• 

• 

• * • • 

• 

• 

• 

• • * • 

• 

• 

• 

• • • f 

• 

0 

0 

0 

I 

-A 
L m 

-A . - 

m-1 

A m-2 ‘ * * 

-A 1 

J 


will have eigenvalues X^ that are the same as the latent roots of A(X) that 
is det[A(X)) = det[A(X)]. Furthermore, it can be shown that the latent 
vectors are subvectors of the eigenvectors of A. If y^ is an eigenvector 
of A for an eigenvalue X^ then y_ L is a subvector of y^. The eigenvector 


of A is given by 
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(3.7) 




V 


, A 
L n» 


V 


0 

-I 

0 


A m-1 A m-2 


lit 
• • t 


0 

0 

f 

* 

-I 

X i I+A l 



J 

'K > 
H- 

1 


Vi 

• 

♦ 


• 

• 

■> Dfr-1* 

<J 

. A i J 


when A^ Is a distinct eigenvalue or A Is not defective. It follows that the 
first n elements of y^ is a latent vector of A(A) . Similarly if z^ is a 
left eigenvector of A then the left latent vector of A(A) will be the last 
..-elements o£ the row vector ^ nnder the seme restriction on \ r 

There is a second relationship between A (A) and A(A) that will be useful 
in the development that follows. The inverse of A(A) in companion form is 
given by 

-1 

- [A(A)]' 1 B adj [A(A)] 


where Badj[A(A)] denotes the block adjoint of A (A) . The block adjoint is 
defined as the adjoint matrix of A(A) with each block matrix, of A (A) treated 
as a scalar element. As an example the block adjoint of A(A) with m - 3 is 


(3.8) 


[A (A) ] 


-1 


AX -I 
0 AI 


0 

-X 


0 

0 


m 


A m-1 A m-2 


-I 

Al+A 


1J 


r» 

-1 

— 

0 


a 2 i+a a-i-a,, 

1 ■** 

A 

a’ 2 i+aa a 

“Aa 2 -a 3 

•> 

i 

i — 
> c 

U5 

AI 

A 2 

-1 

n«i_ 


1 

I 1 

C5 

AI. 

A 2 ! 


0.9) 


Bad j | A (A) J s Bad) 
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The last coltum of B a( jjLA(A)] will always have the form shown with the last 
block equal to 

Consider the inverse of [A(A)) and let A(A) have distinct latent roots. 
The partial fraction expansion of [A(A)] was given by 

i P <n 

(3.10) [A(A)f - l 

i-1 A l 


where P^ Q - (X— [A(X)]~^ evaluated at A ■ X^. It follows from (3.9) and 
(3.10) that the eigenprojector and the latent projector P iQ are related 
since 


(3.11) 


p io - p i0 


• • • 


Pi 


10 

c\ 

.... X i l> 10 


I * • 


x m-l£ 

• * * * A i r i0 


where the first m-1 columns of (3.11) are not important to the development. 

mn 


It was shown earlier that £ P^q «■ I thus it follows from (3.11) that the 


i-1 


latent projectors have the properties that 


mn 

(3.12a) l P = 0 

i=l 

with 


(3.12b) 

mn . 

Jx * p » • 0 

j = 1,2, . . . ,m-2 

(3.12c) 

i *r lp io - 1 



28 


A 

where are the primary latent projectors. If P^ Q is a matrix residue of 
[A(X))’”'*‘ then P^q Is a matrix of [A(X)]"*^. 

The partial fraction expansion of (A(X)] for distinct latent roots 
can be obtained from (3.10) and (3.11) and is 

A 

1 01:1 P 40 

(3.i3) [A(X)f - I rr 

i-i A_A i 


where 


(3.14) 


10 


(X-X ± ) [A(X) ] 


-1 


X-X. 


A 

with P i( j being primary latent projectors. 

It was shown in Section 2 that the eigenprojectors for the A matrix with 
repeated eigenvalues are given by 


(3.15) 


1*^-3 


( V r 1+J+ «! dX . 


{(X-X^ 1 [A(X 1 )f 1 } 


with j = 0,l,2,...,r^ and (•)! = 1 for ( * )j^.O . 

Using (3. 11) , it follows that 

A 

Equation (3.16) agrees with the usual partial fraction expansion formula 
provided that [A(X)] ^ Is a minimum polynomial, factors common to adj[A(X)] 
and det[A(X)] have been cancelled. 

The numerical procedure for computing the latent projectors by the 
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residues will be given for the matrix polynomial 



‘l 0 ' 


r~5 ii 

r 6 -si 

A(X) - 


x 2 + 

X+ 


0 l 


L 1 - 5 J 

L-3 6J 


which has latent roots X^ ■ 1 , ^2 " ^3 " ^4 “ 3 and ^2 * *- nverse °£ 

A(X) Is 

-1 

X-2 


[A(X) l’ 


1 1 

X 2 -5X+6 

-X+3 

1 

(X-l)(X-3) 3 

-X+3 

X 2 -5X+6_ 

(X-l) (X- 3) 2 l 


X-2 

-1 


thus 



(X-l) [A(X)]" 1 



Since X ■ 3 appears twice in the minimum form of [A(X) ] ^ and Si * 1, 
there will be one primary and one secondary latent projector. The latent 
projectors are 



The elgenprojector of A were constructed from the eigenvectors of A in 
Section 2, Since the eigenvectors of A have the latent vectors as subvectors, 
the latent projectors of A(X) can be constructed from the latent vectors. 
Consider the distinct eigenvalues X^ for which the eigenprojectors are 


given by 
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(3.17) 


> -iLfi 

10 z i y l 


Let the order of the matrix polynomial m be 2 then satisfies the equation 


(3.18) 


I -I 
L A 2 X i I+A i 


X £ I -I 


n 

~ A " 


y i 


A 


Vi. 


The left eigenvector z^ must satisfy 




l 

H 

1 

H 
<< 
t 


\ x -i 

(3.19) 

z ± 


■ z i^“^i 

1 



A 2 X^I+A^ 

a 2 

MW ^ 


from which it follows that 
(3.20) z ± * I] 

thus the numerator of (3.17) is 


(3.21) 


y i Z 1 = 


I 


y l Z i^i I+A l 1 ^ 


y i 2 i a i I+ V 


~ y i Z i A 2 


1 i 


x i y i 2 i. 


The outer product z^ y^ is 


(3.22) 


Z i y i " z i^i I+A l 13 


L i J 


y i “ Z i 


dA(X i ) 
“ dX 


y i 


which gives for the distinct eigenvalue or latent root the eigenprojector 
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(3.23) 


10 


dA(X t ) 


i dX y l L 


y i Z l (X l I+ V 


A A 


’Vi A 2 


A A 

y l z l 


A A 

Vi z iJ 


The latent projector, as given In (3.11) Is the (1,2) block of P^q or 

* A 

S y i Z i 


(3.24) 


10 


A dMXj) 


‘l dX 


for the distinct latent root. The elgenprojector for the repeated eigenvalue 
nondefective case was given as 


(3.25) 


10 


j j 

r. yi zi 
yl 7 1 1 

3-1 4 y l 


with the obvious extension to the latent projectors as given In (3.26) 


(3.26) 


10 


V 


y i Z i 


j"l P ^V-l) Aj 

Z i dX y l 


Consider the matrix polynomial 


A(X) 


’ 1 0 ' 


-4.5 0.5" 


"4.5 -1.5" 


x 2 + 


X+ 


0 1 


0.5 -4.5 

— p 


-1.5 4.5_ 


for which the latent projectors are to be found. The latent vectors for the 
latent roots X^ * 1, % and X^ * X^ * 3 are 


A 

» 

H* 

A 

r 

H 

'-l 

1 

H 
L ' 

A 2 

■» ■ 

1 

y l = 

HI 

I 

y 2 = 

1 

y 3 = 

* 1 

1-* 

1 

y 3 “ 

- _1 . 


f 


z l “ [-1 - 1 ) z 2 « [1 - l ] z 2 - [-1 - 1 ] 


- tl -1] 


The two distinct latent projectors are 


A 

p m. 
10 


: Mil 

'1 dX 


>1 -1 


-1 -1 


y 2 Z 2 

4 2 dX y 2 


-2 -2 


-2 -2 


The latent projector for the repeated latent root is 


y 3 z 3 


y 3 *3 


-i MU c 1 P 2 M3J. * 

*3 dX y 3 z 3 dX y 3 


3 -1 


-1 3 


The partial traction expansion of [A(X)] is then 


. -1 -1 -2 2 . 3-1 

IA(X) 1 ' 40,-1) _ x + A (X-2) 2 _ 2 + W-3) 3 


” ith f 10 +? 20 +? 30 “ °> and X lW20 +X 3 f 30 " 1 as re<|uired - 

The formulation of the latent projectors for the repeated latent root 

polynomial when A Is defective In terms of the latent vectors remains as a 

problem. Attempts to formulate the latent projectors for the defective 


ease have been unsuccessful. Future work will be devoted to this problem. 





4. Projectors and the Sign Matrix 


It was shown In Section 2 that the primary elgenprojectora for the 
mnxmn A matrix are given by 


(4.1) 




1 — 


where is a diagonal matrix with ones along the diagonal of the ith 
pseudo-* Jordan block. The secondary eigenprojectors were defined as 


(4.2) P Aj - Q F i;J Q" 1 j - 1,2 


when A has generalized eigenvectors for the eigenvalue 

Assume that A has q^ elgenvalueswlLh Re(A^)>0 and q ^ eigenvalues having 
Re(A i )<0 and no eigenvalues with Re(A^) * 0 so that q 1 +q 2 “ q* bet P + be 
denoted as the sum of the eigenprojectors with Re(A i )>0 and P“ the sum of 
the eigenprojectors with Re(A i )<0; that is 


(4.3) 


(4.4) 


P* * 


1-1 


10 


Q l 1 F 


1-1 


10 



p" 


q 

l 

i-q^+l 



Q l 

i-qj+l 




where It has been assumed that the first q^ eigenvalues have Re(A^)>0. 
The sign of a matrix, denoted by S, will be defined as the matrix 


(4.5) S « Q{sign[Re(A]}Q“ 1 
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where Re (A) denotes the reel part of the eigenvalues of A or the diagonal 
elements of the Jordan matrix J, Let E^ be the ranXran matrix with diagonal 
elements equal to 1 If Re(X jL )>0 and zero for Re(X i )<0 and E 2 be the com- 
plement to such that Ej+E,, - I. The sign matrix can then be defined as 


(4.6) 


s ■ Q E Q’ 1 - Q E, Q" 1 


Q[ t F 10 


1-1 


? , -1 

V F 10 3Q 
i-q x +l AU 


therefore S Is equal to 


(4.7) 


S - P + - P' 


Knowledge of the elgenprojectors is sufficient to construct the sign matrix, 
similarly it can be shown that knowledge of the sign matrix is sufficient 
to construct P and P . Assume that P is given as 


(4.8) p+ - \ (S+I) - | [QE 1 q" 1 -QE 2 Q“ 1 +QQ" 1 ] 


but since ■ I then 


1 (4.9) 


P + - \ Q[E 1 -E 2 -fE 1 +E 2 ]Q‘ 1 


QE i Q ’ 1 ' q Ji W 


-1 


It Is not difficult to show that I’ is given by 


( 4 . 10 ) 


P” = | U-S] 
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The computation of the sign of A is a rather simple task, Roberts [10], 

gave an iterative algorithm to compute S which is based on Newton's method 

2 

for computing the square root of S * I, The algorithm is 

(4.11) S(l+1) - \ {SUHISU)]” 1 } S(0) - A 

where the index 1 denotes the 1th iteration. The algorithm will converge 
quadratlcally to S provided that A has no eigenvalues on the jw axis. The 
simplest test of convergence of (4.11) to the sign of A is to compute the 
trace of S 2 (l) at each iteration. Since S 2 (i) converges to X, then trace 
[S] will he sm. 

Several accelerated versions of (4.11) have been described in the 
literature, Roberts [10], Hoskins and Walton, [11], and Mat they s, [12]. 
Numerous applications of the sign algorithm to system analysis have been 
given in the literature, [13]- [16], 

The example below gives the sign of A where A is 

-1 -1 9 -3 

-3 1-1 7 

9 -3 -1 -1 
-1 7-3 1 _ 

with eigenvalues X^ “ 1, m 3, X^ " X^ « -2 and ^ ■ L The sign of A is 

0 0 I 0 

„ 0 0 0 I 

!•> 

1 0 0 0 

0 I 0 0 
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which converged in 5 iterations. The trace of [S(i)] 2 was 


iteration trace 


1 6.90278 

2 4.3857 

3 4.01697 

4 4.00006 

5 4.000P0 


The positive and negative projectors, 


and P , were found to be 


P 


+ 



1 0 
0 1 
1 0 
. o 1 


1 0 
0 1 
1 0 
0 1 


m P 4- P 
10 


20 



The sign matrix of A when all eigenvalues have Re(A^)>0 will be I where- 
as the sign of A with Re(l^)<0 will be -I. Eigenvalues along the jto axis 
can be removed from the axis by an origin shift or by computing the sign of 
(A+pI) where p is a real number. The eigenvalues of A+pI will be X^+p since 
A+pI = Q(J+pI)Q \ Eigenvalues belonging to a Jordan block cannot be separated 
nor can eigenvalues along the ju) axis be split by the described procedure. 

A method of separating eigenvalues according to their magnitude is to 
compute a new matrix by the bilinear transformation 


1 
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(4.12) 


A q - (A-pI) (A+pi) 


-1 


where p has the same meaning as above. All of the eigenvalues with | X ± | <p 
will be mapped Into the half plane Re(X i )<0 with others mapped Into the plane 
with Re(X i )>0. This procedure is more general than the origin shifting 
method since the spectrum splitting will be according to the magnitudes of 
^i* ^ ie *' wo •“skhods, shifting and splitting can be combined if desired to 
isolate any circular region of the eigenvalue space. For example, the matrix 
Aq given by 


(4.13) 


A q - (A+Pjl-Pgl) (A+p x l+p 2 l) 


-1 


can be used to isolate eigenvalues inside a circle of radius p 2 centered at 


Pi 


1* 


As an example of the bilinear transformation procedure, let 


3 

>2 

-9 

6 


-2 -9 
3 6 

6 3 

-9 -2 


6 

-9 

-2 

3 


trace - 6 


■ i 

4 


which has eigenvalues X^ * -1, X 2 “ 2, X^ = -5 and X^ « 10. If the value 
of p “ 4 is selected then all eigenvalues inside of the circle p » 4 will 
be mapped to the left half plane and those outside the circle will be in 
the right half plane. The sign of A is 


S K 


0-10 0 
-10 0 0 

0 0 0 -1 
0 0-1 0 
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with projectors 


1 

2 


1 -1 

>1 1 


0 

0 


0 

0 


0 ol 
0 0 ! 

1 -i 
•l i 


i 

i 

0 

0 


1 

1 

0 

0 


0 

0 

1 

1 


0 

0 

1 

0 


The projectors 1? and P are idempotent matrices and satisfy the proper- 
ties of the eigenprojector given earlier. It is not difficult to show that 
the positive and negative projectors can be used for the spectral decomposi- 
tion of A. Since P* -QE 1 (f 1 with E x - diag[l 111. .000] with die 
ones in the first locations. The product AP will have the eigenvalues 
with | | > p whereas AP will have eigenvalues | | <p with all other eigen- 
values zero. Using the example 



5 

5 

-15 

15 




1+ l„+ 1 

A 3 AP 

4 

-5 

5 

15 

-15 

trace 

* 5 “ 

A 4 -A3 

1 -15 

15 

5 

-5 






_ 15 

-15 

-5 

5 

l 





* 

1 

1 ■ 

-3 

3 





A" - AP~ - ■“ 

1 

1 - 

-3 -: 

3 


trace 

* 1 * 

^2~^1 


-3 

-3 

i i 






3 

-3 

i i 
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i, ^ r 

The sura of AP and AP must be A since P^+P - « I. 

HH - 

It is obvious that P and P can be decomposed into eigenprojectors for 
the eigenvalues A^, X2» ^3 ant * A^. The details for the determination of the 
eigenprojectors have been covered in Section 2 and will not be covered at 
this point. 

4- M 

The positive and negative projectors, P and P , have been defined in 
(4.3) and (4.4). It follows from (4.3) and (3.11) that P^ is also given by 


(4.13) P + « j 1 P 
i=l 


iO 


j£ p io<V+V 


n 1 

T 


I 1 ? 


i»l 


10 


l i=l 


10 2 




1=1 


i 10 


for the companion form when A(A) is a second-order polynomial. 
Slmilarlly, P = I-P + thus 


(4.14) P" = -l 


,i=q ^+1 


10 


-l 

i-q^l 


i=qj+l 


P iO (X l I+ V 


P :L0 A 2 


-l 


i-q^+1 


iO 


-l 


i-q +1 


A.P.. 
x iO 


The secondary eigenprojectors and latent projectors are not needed in the 
decomposition of A into A and A . This can be shown from the definition 
of A + and A . 

The individual eigenprojectors P^q can be computed by repeated use of 
the sign algorithm. Assume that n = 4 with = 1, A£ = 2, A^ = 3 and 
A^ = 4. The bilinear transformation can be carried out first with p = 1.5 




AO 


to separate from the other eigenvalues. If denotes the sign of A Q1 
with 


(A. 14) A 01 - (A-1.5IKA+1.5I)” 1 


then will be given by 


(A. 15) 


P 10 ■ f <I - S 1 ) - ''i 


I 

with P - £ P iQ . The next step in the procedure is to compute A n0 with 


i-2 


02 


P 2 “ 2,5 and compute the sign of A Q2 . The negative projector of Aq 2 will 


be the sum of P^q and P 2 q or 


P 2 “ P 10 + P 20 “ 2 • 


Thus P 2Q = l> 2 “ P io' Tills process can be continued until each eigenprojector 
lias been found. Since the eigenprojectors also give the latent projectors 
as the upper right block, the latent projectors will also be known when A 
is in companion form. 

A method of computing the projectors of a matrix has been discussed in 
this section. It has been shown that the eigenprojector for any general 
matrix can be computed from the sign of a matrix. If A is xn companion 
form, the latent projectors can also be found from the eigenprojectors. 


5. Solvents of Matrix Polynomials 


The concept of matrix polynomials was Introduced in Section 1 of this 
report where it was shown that A(X) arises when n simultaneous equations of 
mth-ordar are used to define the time-behavior of a dynamic system. The 
eigenprojectors for a matrix in companion form was discussed in Section 2 
and the latent projectors were described in Section 3. The application of 
lambda matrices to the dynamics of systems has been described by Frazer, 
Duncan and Collar, [ 17 ] # in their book El ementary Matrices and Some Applica- 
tlons to Dynamics and Differential Equations. The concept of solvents or 
matrix roots of a matrix polynomial will be given in tills section. It will 
be shown in the next section that solvents are useful in solving sets of 
differential equations. 

Let A(X) be defined as a mth-order matrix polynomial with nxn matrix 
coefficients. The associated mnxmn A matrix, which will be called the 
block companion matrix, is given by 


(5.1) 


0 I 0 ... 0 

00 I ... 0 


0 0 0 

-A -At - A 0 
m m-1 m-2 



The eigenvector matrix Q for A will always have the form 




" Q 1 

q 2 




Q 1 A 1 

q 2 a 2 

Vm 

(5.2) 

Q = 

2 

2 

2 



<¥i 

q 2 a 2 

Q a 

m m 



_vr l 

• • • ■ 
• i • 

Q 2 A 2 

Q A m "’ 
m m 
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when A has distinct eigenvalues or A(X) has distinct latent roots. The sub- 
matrices Qj will be a matrix of latent vectors y^ for the latent roots 
It will be assumed that exists and is Invertible; under the above assump- 
tions, Rj “ 3 so ^ vent matr *- x polynomial and satisfies the 

equation, [1], 


(5.3) 


R 1 " + A. r " 1 " 1 + 

.1 1 j 


+ A 


m 


j » 1 , 2 , . . . ,m 


The proof of this is straightforward if A(A) is considered. Let A(A) be de- 
fined as the block matrix Al-A; it then follows that 


1 

• • , > 

u. 

1 

• O Hi 

. . o 

. . o 

• • o 

J 


Vj 


’ 0 * 

0 

* 

• • • • • 

0 0 0 A. -I 

j 

A A , A o A„ A.+A, 
m m-1 m-2 2 j 1 


# 

• 

Q.A"" 1 
3 3 


■ 

0 

m * 


where (5.3) is given by the last row of (5.4) provided that A^Q^ = Q^A,. 
must hold from the first row of (5.4). 


The block matrices A. are nxn diagonal matrices constructed from a 

j 

subset of the latent roots of A(X) or the eigenvalues of A. Each R^ is de- 

A 

fined in terms of n latent roots and n latent vectors y^, as an example, 
consider the matrix polynomial A(X) with 


*1 o’ 

o 

- «— 
-5 2 


r 

in 

\ 

r- 
l 


x + 


X+ 


0 1 
L J 


2 -5 


-5 7_ 


A(X) = 


43 


which has latent roots X^ - 1, X 2 ■ 2, X^ ■ 3 and X^ - 4. The right latent 
vectors of A(X) are 



The solvents are constructed from the latent roots and latent vectors by 
forming the 2x2 matrices where * {y^} such that is invertible 

with j indicating a subset of latent vectors. Noting that is 

singular, then 

A A A r 2 -i ' 

R 1 = ty l y 3^ diaglX l X 3^ y ] y 3 1_ = i 2 

3 _i 

R 2 = ty 2 y A ldiag[X 2 X 4 ]ty 2 y 4 1_1 = 3 

It can be shown that and R^ satisfy (5.3). 

The eigenprojectors for the companion matrix A are given by 


(5.5) 




• 


• 

• 

I 

• 


• 

• 

x.i 





i 

• 


• 

• 

, 2,. 
XjI 





i 

• 


• 

• 

• 

• 


• 

• 

• 

• 


• 

• 

X™ _1 I 

1 <N 


where the first m-1 block columns are not important in the development. 
Suppose that the latent projectors are known and A has distinct eigenvalues. 
It can be shown that if is defined as 
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(J+l) n A 

(5.6) 

Q J +1 ' i"Jn+l Pl ° 

and 



(j+Dn 

(5.7) 

Q4+1 - 1 . 

3 3 i-jn+l 


j ** 0|lj2| i • • |in— 1 


*10 


then Rj Is defined by the latent projectors of A(A) provided that Is in- 
vertible. If Qj Is singular, the latent projectors are reordered until a 
set is found for Q"'*' to exist. 

j 

The solvents for repeated roots can be defined by noting that A. may 

J 

Include several Jordan blocks. Assuming that the multiplicity r^ is less 
than n, and that full Jordan blocks are included in A., Q. Is defined 

j j 

as in (5.6). The product QjAj must be modified since Aj is no longer 

diagonal but may include the ones on the super diagonal of A. due to the 

j 

included Jordan block. Assuming that the Jordan block is defective, then 


QjJj will be defined as 


(5.8) 


C+ 1 )" . 

V.1 ■ 1 IX i P iO + P il I 

3 3 i=jn+l 1 1U 


The matrix polynomial 


A(A) 


"1 o' 

0 

"-3.5 0.5' 


2.5 -0.5" 


x> 


A+ 


0 1 
4 -* r 


_ 1.5 -4.5_ 


3 ,5 5.5. 


has latent roots A^ = 1, <= 3, A 3 

block for A = 2 is then 

2 l' 


= 2 and A^ = 2 with = 1. The Jordan 


J 3 = 


0 2 
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The latent vectors of A(X) are 


A 

m ”1 

1 

A 

1* 

*1 

■ ■ 

1 

A 2 

1 ' 

y l " 

1 

w • 

y 2 " 

-1 

y 3 " 

1 

► a 

y 3 

“I. 


^ * [1 0] z 2 » [1 -1] zjj - [1 -1] z\ - [-1 -1) 


The first solvent Is given by 


R i “ WI 1 


„ A £ A _1 

^ P ioV,L P 10 J 


1-1 


1-1 


2 -1 
-1 2 


which gives A(R^) = 0 as required. The second solvent requires the primary 
latent projector and the secondary latent projector which are 


with 


A 1 

1 1 

n 1 

-1 1 

P m —• 

*30 2 

3 -1 

a. m 

’ Q 2 P 31 “ 2 

-1 1 
• m 


Qo J 


2 2 


P 30 X 3 * P 31 or 


Q 2 J 2 


1 

2 


-1 

-1 


1 

*T 


1 3 
5 -1 


The second solvent is 


„ T „-l 

1 3 


H* 

1 

-1 

1 

r 5 -ii 

R 2 “ Q 2 J 2 Q 2 

5 -1_ 


I 

Ui 

i 

H-* 

1 

" 2 

cn ’ 
I 


which gives A(R 2 ) = 0 as required. 

The extension to higher order polynomials Is a simple matter and does 
not require additional analysis. Each solvent, R^,R 2 » * * * > R m > * s ^ ound by 
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the procedure given In the preceeding work. 

The solvents can be determined by the sign algorithm provided that the 
spectrum of A has the required distribution. The first step in the proce- 
dure Is to establish A for the matrix polynomial and assume that the 
multiplicity of the repeated eigenvalues are less than n. Let AjJ , . . . f | A | 
have magnitudes less than p^, and compute the sign of A^ where 
Aj)i « (A~pjl) (A+Pjl) . The sign of A Q ^ will have n eigenvalues of -1 and 
mn-n +1 eigenvalues. The sign of Aj^ can be arranged In the form 

(5.9) « Sign <A D1 ) ■ Q 


-1 


n*n 


vl 


mn-n*mn-n 


by row-column Interchanges of the sign as computed by (4.11). Let the 
eigenvalue matrix (5.9) be denoted by then 


(5.10) (S i +J H ) “ [QJ I1 +J I1 +J 11 Q ^ 


-1 


Suppose that Q Is partitioned as 


(5.11) 


Q - 


*1 Q 12 


Q 21 Q 22. 


Q 1 A 1 


Vi " 1 


Q 2 A 2 


vr 1 


• i • 

• ■ • 
» • • 


m 

Q A 

m ra 


Q A™ -1 
^m m 


which when substituted into (5.10) gives 


ts 1+ Jn3 - 2 


Ql o 


0 Q 


22 J 


.-1 


(5.12) 
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The similarity transformation then gives 


(5,13) 




«1 ° 


0 Q, 


A 1 0 


0 L 


q” 1 0 


0 <L 


R 1 0 


0 R„ 


The similarity transformation T required to reduce A to the block diagonal 
form is given by 


(5,17) 


^^nrM, R -(mKL) 

, . ■ v urH_-(mfl)+l 
(-1) R m 

^_ 1 )iirfl R -(mfl)+2 

m 


„m-l D m-2 
R 1 R 2 


(-1) m ^ 1 l 


where R^ is a solvent of A(X). 

The spectral decomposition of A will now be shown using the second ex- 
ample in this section. Let A(A) be defined as 


0 1 


1 0 0 -3.5 0.5 

X l + A+ 


1.5 -4.5 


2.5 -0.5 
-3.5 5.5 


with A^ = 1, X 2 ® 3, Ag = A^ ** 2 and ^ « 1. The companion form A Is 
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0 

0 

1 

*» 

0 

A - 

0 

0 

0 

1 


-2.5 

0.5 

3.5 

-0.5 


_ 3.5 

-5.5 

-1, 5 

4.5 _ 


with solvents 


R, 


From (5.17), T Is 


2 -1 

-1 2 


R„ 


2.5 -0.5 
0.5 1.5 


with 


T 


1 0 

0 1 

2 -1 

-1 2 


3/8 

> 1/8 

■1 

0 


1/8 

5/8 

0 

-1 


A m “ A„ 

m l 


T A T” 1 


2-10 0 

-12 0 0 

0 0 2.5 -0.5 

0 0 0.5 1.5 J 


* 

t 


i 


2 

The characteristic equation for the upper block is X -4X+3 and the lower 
block has X^-4X+4 as its characteristic equation. The eigenvalues X^ * 1 
and X 2 = 3 are in the upper block with X^ * 2 and X^ ** 2 in the lower block. 
The computations for this example were checked by the sign algorithm with 
a shift of -1.8 and p x - 0.5. The two solvents, ^ and R 2 computed by the 
sign algorithm, agree with the values given. 



6. Solution of a System of Differential Equations 

The mathematical tools developed in the previous sections will now be 
applied to the time domain analysis of systems. Assume that the system has 
been characterized In the first-order form, usually called the state vari- 
able form, with states x(t) such that x(t) satisfies the differential equa- 
tion 

(6.1) iiy . p) - Ax(t) Bu(t) 

where A is a general nnXran matrix, B is mnxkn, x(t) is a maxi vector and 
u(t) is knxi. Let Z(t) denote the outvector with 

(6.2) Z(t) - Cx(t) 

where Z(t) is kmxl and C Is kmxmn. The vector u(t) will be considered as 
the input to the system or a control vector if the system Is a control sys- 
tem. It will be assumed that A,B and G are constant matrices which will be 
referred to as the system triplet. 

It will be assumed that the system is stable, all eigenvalues of A 
have Re(Xj)<0 except for distinct eigenvalues along the jw axis and multiple 
eigenvalues at the origin, X^ * 0. 

The solution to (6.1) car* be expressed as 

(6.3) x(t) * 0(t,t Q )x(t o ) + f 0(t,T)Bu(T)dt 

*0 


49 


50 


where 0(t,t Q ) denotes* the state transition matrix, (STM). The state- 
transition matrix satisfies the differential equation 

d0(t,t ) 

(6.4) - — - A 0(t,t Q ) 0(t o ,t o ) - I - 0(t,t) 

where Q(t,tg) Is ranXmn and I is the Identity matrix. There are numerous 
methods of obtaining 0(t,tg), several of the methods will be described. 

Since A Is a constant matrix, the Laplace transform of (6.4) can be 
taken with 


(6.5) L[~* 0(t,O)]« s0(s)-0(O,O) - s0(s)-I 
thus the transform of (6.4) Is 

(6.6) 0(s) = [sI-A]" 1 - QtsI-Jj’V 1 


where J is the Jordan form. Noting that s in the Laplace domain is equiva- 

— 1 " 5 1 

lent to X in the eigenvalue domain, 0(s) is equivalent to [XI-A] = [A (X) ] . 

It therefore follows that 0(s) can be expressed as 


3 p -n K 

(6.7) 0(a)- l h£-+ J 1 

1-1 s S 1 J-l 

The inverse Laplace transform of 0(s) 


— i*-r] 

(s-Si ) 3 

is L 1 [0(s) ] or 


( 6 . 8 ) 


0(t,O) = Q{exp(Jt)}Q 


-1 


q 

I 

i=l 


{P iQ exp(s jL t) + 


l. 


d=i 


p. .t J 

ji 


exp(s^t) } 


*The normal use of $(t,tg) as the state- transition matrix will not be made 
due to the use of (j> as the modal matrix, [17]. 
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where Is a eigenvalue of A or a root of the characteristic equation 
det(sI-A] m 0. 

The Laplace transform is usually taken with t Q ■ 0 rather than on arbi- 
trary value of tg. The state transition matrix 0(t,t Q ) can be found by us- 
ing the semigroup properly 


(6.9) 0<t,t o ) - 0(t, 0)0(0, t Q ) 


,-1 


where 0(0, 1^) - [0(tg,0] if t^>0. It follows that ©(t,!^) is given by 


( 6 . 10 ) 


0(t»t Q ) * Qtexi>[J(u-L 0 ))}(r :L - l IP Q exp(s i (t-t () )] 

v 1=1 




i 


+ X 1 exp [s. ( (t-t n ) ) } 


i-1 




The analytical solution to (6.1) can be expressed as 


(6.11) x(t) » Q exp [J (t-t Q ) )Q“' L x(t 0 )+Q f exp[J(t-T) )Q _:L Bu(x)dX 

+ t* . _ 


where Q is the eigenvector matrix and J is the Jordan form. If the eigen- 
projectors P^q and P^_. are used, then 


( 6 . 12 ) 


II p (t~t )^ 

<(t) = i {P i0 exp[s i (t-t 0 ) ] -I- }’ x — expls^ (t-t 0 ))x(t 0 ) 


q 

i 

i-1 

ft 3 


t Q i-1 


n & d 

l {exp[s < (t-'C)]P <0 + l L exp[s 1 (t-T)P ij }Bu(T)dt 


10 j=l 3* 


where the order of the operation under the integral has been reversed for 
convenience. This reordering is permissible since exp(*) is a scalar. 
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The system defined In (6.1) Is said to be controllable If the Input u(t) 
(or the control) can drive the Initial states x(tg) to the origin x(t) « 0 
(or to an arbitrary value x(t) ** x f ). The usual test of controllability is 
to examine the matrix 


(6.13) 


Q c - [B,AB,A 2 B A 1 ^) 


to determine if Q is invertible. The system is controllable if det Q + 0 
c c 

or if Q c is invertible. Controllability of a system can also be measured by 
considering the product P^B for 1 - l,2,...,q and j » l,...,#.^ for the 
eigenprojectors. The system mode exp(A^t) is not controllable if 


(6.14) 


P, .B ■ 0 
ij 


for the primary and secondary eigenprojectors of the eigenvalue This 

test implies that u(t) cannot drive the mode exp(A^t) as the mode is follow- 
ing the natural response rather than a forced response. 

The partial fraction expansion method of determining 0(t,t^) for all t 
is not a computationally efficient process and would be used only when the 
analytical form of 0(t,tg) is desired. For computational purposes, the 
state-transition matrix 0(t,O) can be determined more efficiently by ex- 
pressing 0(t,O) in the form 

00 j *j 

(6.15) 0(t,0) = exp At = X + J 

J-l j 

The series is then used to compute 0(t,O) for a small value of t with trun- 
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cation of the series when the change In 0(t,O) is beyond the accuracy of the 
digital computer. Assume that the small t Is taken as At, with 0(At,O) de- 
termined from (6.15). The semigroup properties is utilized to find 0(t,O) 
for t « kAt by the operation 

(6.16) 0(t,O) * [0(At,O)] k - 0(At,O) 0(t-At ,0) 
when A is constant. 

The solution vector x(t) depends upon the integral as well as the state 
transition matrix. The expression in (6.3) can be given in a recursive form 
which is more convenient and which has computational advantages. Let (6.3) 
be written as 

(6.17) x(t) = 0(t,t Q ) x(t Q ) + r(t,t Q ) 

where 0(t,tQ) has been computed by means of (6.10) and where T(t,tp) is the 
integral. The solution vector at t+At can be expressed as 

(6.18) x(t+At) = 0(t+A,t)x(t) + r(t+At,t) 

Substituting (6.17) into (6.18), 

(6.19) x(t+At) « 0 ( t+A t , t) 0 (t , t Q ) x( t Q )+r ( t+At , t)+0 ( t+At , t) T ( t , t Q ) 
or 

x(t+At) = 0(t+At,t^)x(tQ) + r(t+At,t Q ) 


( 6 . 20 ) 
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Knowledge of r(t+At,t) is required along with 0(t+At,t) or 0(At,O), since A 
is constant, to find x(t+At). The vector F(t,tg) satisfies the equation 

dr(t,t n ) 

(6.21) ^ - A r(t,t Q ) + Bu(t) r (t 0 ,t 0 ) - 0 

and can be found a 4th-order Runge-Kutta algorithm. Since F(t+At,t) is re- 
quired for (6.20), the value is found from (6.21) by integrating (6.21) from 
t to t+At for all t. Equation (6.19) is then used in a recursive manner 
with 0(t+At,t) *» 0(At,O) for constant A. The series given in (6.15) is 
used to compute 0(At,O) which remains constant thereafter. 

If u(t) is slowly varying and At is small, a reasonably good approxi- 
mation may be obtained from the z-transform, [3], [18], of (6.1). The 
Laplace transform of (6.1) is 

(6.22) [s'I-A]X(s) = Bu(s) 

where the initial condition x(0) has been neglected. Equation (6.22) can 
be written as 

(6.23) x(s) = [sI-A] _1 Bu(s) = G(s) U(s) 

Taking the z-transform of G(s) witli a zero-order-hold gives 

_ ST 

(6.24) G (z) = Z(i— — [sI-A] _1 ]B 


which is found to be 


(6.25) 


G(z) » (zI-exp(AT) ]”^j.exp (AT)-I]A -1 B 
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where T - At. The associated difference equation of (6.1) is obtained from 

(6.25) with* 


(6.26) x(k+l) * [exp(AT)-I]A ^B U(k) + exp(AT)x(k) 

with t * kl’, (k+l)T “ t+At and x(0) the initial condition vector. The 
validity of (6.26) can be shown by noting that 


(6.27) 0(t,O) - exp (A t) 


(6.28) 


r T fT 

0(t,T)B dl = exp[A(t-T)]B dT 
•'0 •'O 

= [exp(AT)-I]A ^B 


Equation (6.26) is valid provided that the assumption of slowly varying 
u(t) and small At is not too strict. 

The eigenprojectors can also be utilized to separate the solutions into 
modes if desired. Recalling that [sI-A] ^ is equivalent to [Al-A] ^ and 
that [AI-A] ^ can be expanded as a partial fraction, (6.24) can be refor- 
mulated as 


(6.29) 


G(z) = [A(O )]" 1 + l (■ 


z— 1 


P i° + l± 


P T J 

_±i 


• 1 s, 1 SjT , SjT.j+l 

;l=1 i z-e 1 j=l (z-e y 


J8 i\ 
e ) 


which can be written as a series of transfer functions G, . with 

ij 


The argument of x(') and u(-) is (k+l)T with T being dropped for convenience 
hereafter. 






1 


I 

i 


I 

i 
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[g oo + yV oo] 

0 j-i 13 


q A 

+ l [G Q (z) + i G. . (z) ]B U(z) 
i“l 1U j-1 


» Equation (6.31) defines a set of transfer functions that can be placed in 

parallel with each having u(z) as an input and each block can be implemented 
i 

j separately by a matrix difference equation. 

Recalling that the eigenprojectors are defined by the right and left 

i 

i- 

eigenvectors , the eigenprojectors need not be stored in the digital computer. 
[ The eigenprojectors are stored and used in (6.29) to construct (6.29) at 


each T. 

The discussion in this section has been based on A having a general 
form. Assume that A is the companion form of the matrix polynomial A(A) 


with 


! 

| | (6.30) G(z) - G Q0 (z) + 

' The vector x(z) is then given by 

Pi I 

r 

(6.31) x(z) - Gqq(z)B U(z) 


,m A ,mrl A 

(6.3 2) I-2+A 1 - J ^ + ... + A m x= u(t) 


The canonical form of x(t) is the same as given in (6.1) but with 


(6.33) i(t) = 

?t(t) 

x(t) 

• 


OX 0 ... 0 

00 X ... 0 

• • • • • • 


x(t) 

A 

x(t) 

• 

1 

i 

• 

- .9 


• . • • • • • 

* • • ♦ • • # 

-A -A * -A ,y • • • -A . 

m nri m-/. .1 


• 

• 

• 
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+ . u(t) ■ A x(t) + B u(t) 


Rather than use (6.33), the solution vector x(t) can be determined directly 
from (6.32). Taking the Laplace transform of (6.32) gives 


(6.34) 


[Is m +A.. s m ~’*' + ... + A .jtf-A ]x(s) ■ U(s) 
1 m -1 m 


where the initial conditions have been neglected. Equation (6.34) can be 
rewritten as 


(6.35) X(s) » [A(s)f U(s) » G(s) U(s) 


which can be expanded into a partial fraction. 

The z-transform of (6.35) can now be found provided that the zero-order- 
hold is included and u(t) is slowly varying. Letting 


(6.36) 


-sT 

G(s) - G(s) 


gives 


(6.37) 


— c?T — sT 

G(z) = 5S O 3 —-- - [A(s) f 1 } “ G(s) } 

b b 


The z-transform of (6.37) is 
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(6.38) 


q a 

G(z) - [A(0) J” 1 + l {- ^ 


i-1 8 i 


8 i T 1-1 
(z-e 1 ) J 1 


A 4 

L P 44 T3 e 1 

+ l - iL T-S ) 


(x-.“ lT ) J+l 


from which the difference equation for x(k+l) can be determined. 

The expression in (6.38) would be Implemented on a digital machine in 
parallel form rather than as a series representation. As an example , let 



with roots A 1 « -1, A 2 ■ -2, A^ ■ *-3 and A^ » -4. The latent projectors are 


• 

1 

Pi 

1 

A 

P _ * 

-1 - l ’ 

A 1 

P = i 

1 -1 ' 

_1 

1 _ 

20 

-1 -1 

> , * 

30 2 

-1 1 

w • 


/ 



with A(0) obtained from A(s). The transfer function G(z) is 


where 


G(z) - G 00 (z)+G 10 (z)+G 20 (z)+G 30 (z)+G^ 0 (z) 


G 00 (z) " 24 


7 5 
5 7 


G 10 (Z > = 


z— 1 

' 1 1 

1 

’ 1 l" 

2(z-e" T ) 

_1 1_ 

“ 2 

1 1 
■*. -J 


1-z 


-1 


1 -T -1 
1-e z 


G 20 (z) ~ ' -21 r ‘ 
2 (z-e ) 


r-L -1' 

1 

-1 -1 ‘ 

t — 
1 

i-* 

1 

t— 

■ 

a 2 

_-l -1 _ 


1 - 1 

1-z 

. -2T -1 

1-e z 
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G 30 (Z) 


z-1 


r i -ii 


2(z-e“ 3T ) -1 1 


1 -1 

-1 1 


l-z 


-X 


, -3T-1 

1-e z 


V z > ■ TTT?^ 


-1 i 

2<z-e~’ ,i ) L 1 -1 


1 

2 


l-z 


-1 


-1 l] 

1 -lj l-e- 4 V l 


The implementation as a difference equation would then be 


with 


X 0 (k+1) ■ 24 


7 5 


1 

f u(k+l) 


L5 7J 


X^k+1) * 


1 1 
1 1 


(u(k)-u(k-l)) + e” T X^k) 


X 2 (k+1) - \ 


•1 -1 
■1 -1 


(u(k)-u(k-l)) + e“ 2i X 2 (k) 


X 3 (k+1) - ^ 


1 -1 


-1 1 


— w A 

(u(k)-u(k-l) + e X^(k) 


X 4 (k+1) - j 


-1 1 


1 -1 


(u(k)-u(k-l)) + e" 4T X^(k) 


X(k+1) = l X, (k+1) 
k-0 K 


The parameter T is the integration step size (or sampling rate) with t = kT. 
Different step sizes can be selected for the modes cxp(-s/l‘) if better 
accuracy is desired for the higher frequency modes. 
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The use of the eigenprojectors and the latent projector for solving the 
state equation of (6.1) has been described. The extension to the solution 
of n rath-order differential equations has been given. 


7. Spectral Decomposition of Differential Equations 


The increasing complexity of modem systems generally requires a large 
number of algebraic operations on system equations and the corresponding solutions to 
these uquatlons to characterize the dynamics of the systems. It is there- 
fore essential that the analysis of large scale systems be carried out on 
subsystems or subsets of the equations. A dissertation by K. S. Yoo, [15], 
used the concept of mode decoupling to analyze an optimal control system. 

The overhead for the decoupling procedure was moderately high but the total 
computational task for the analysis was decreased when compared to the 
usual procedure. Popeeva and Lupas, tl6], published several papers describ- 
ing the procedure, the referenced paper was the first reference to Yoo's 
work In the open literature. The original idea for decoupling, or order 
reduction, was probably due to Roberts, [10]. 

Consider the system equation given in Section 3 with 

(7.1) , Ax(t) Bu(t) 

Let T denote the mnxmn transformation matrix 

(-l)V mfl 

m 

... miV 5 

m 

H) V»« 
ra 

• * • • 

♦ # ♦ • 

... X 


where I is the nXn identity matrix and is the solvent of 


(7.2) 


T = 


( 2 ) 


m-1 


m-1 


-R 

-I 

R„ 


-1 


, ro -2 


R 


-2 


R' 


w3 
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f 


(7,3) 


R® + A.R?" 1 + ... + A ~R. + A - 0 

i 11 wl 1 D 


Define the vector v(t) by the transformation 


(7.4) 


v(t) * T~ l x(t) 


or x(c) » X v(t) . Assume that A is the companion form for the mth-order 
differential equation, It follows that v(t) must satisfy 


(7.5) - T" 1 * T v(t) + a "’ 1 B u(t) 


Consider the matrix t" a AT when m ■ 3 and T Is as given In (7.2). It follows 


(7.6) 


A b - T -1 AT 


a bi 0 0 

0 A B2 0 

0 0 A„, 


*1 0 0 
0 R 2 0 

0 0 R„ 


or in general for the mth-order differential equation 


(7.7) 


Ag - diag[R 1 ,R 2 ,R 3 ,...,R m i 


Recalling that the solvents R^ were given by 


(7.8) 


- _i (i+Dn a A ( j+1) n A , 

\ ■ Vk\ ■ 1 j Jn+1 lX l P iO +P il H J in+1 V 


Ag aitd T are completely known provided that the latent projectors are known. 


I 


Tha solution to (7,5) Is obtained from the m equations 
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dv. (t) 

(7.9) ■ R ± v(t) + B iU (t) 

where 


(7.10) 


B » T^B - 


i. m J 


with Initial conditions v^(t Q ) obtained from (7.4) with x(tg) given. The 
mnxi vector x(t) as given In (6.33) is 


(7.11) 


x(t) “ 


x(t) 

x(t) 


a 



and x(t) ■ T v(t), The solution to the rath-order differential equation of 
(6.32) is 


( 7 . 12 ) 


x(t) 


( 2 ) 


m-1 


(v 1 (t)-R2 1 v 2 (t) + R^ 2 v 3 (t) - ...} 


The decoupling procedure described is useful in solving a large set of 
differential equations as it allows the set to be decoupled into a number 
of equations, A mth-order differential equations with nxn coefficients can be 
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solved by generating ra first-order equations with n><n coefficients. These 
equations can be reduced further if desired by decoupling any of the m 
first-order to several first order differential equations with coefficients 
less than nxn. The limit to the decoupling procedure will be mn first-order 
equations. 

As an example of the above let A be the general matrix 


-168 

48 

279 

-218 

8 

50 

-136 

32 

233 

-180 

14 

36 

-80 

16 

95 

-58 

8 

18 

-60 

12 

75 

-50 

8 

14 

-40 

8 

50 

-32 

3 

10 

-20 

4 

25 

-16 

2 

4 

;*1 = 

1 » *2 

- - 2 , 

X 3 “ 

-5. X 4 

C5 - 


Ag « -40. This matrix can be reduced to the companion form by the Krylov 


transformation 


(7.13) A c = K A K _1 


where the structure of K is given in the Appendix. The associated matrix poly- 
nomial A(A) is 



i o " 


’ 702.051 

- 734 . 624 " 


7688.79 

- 6761.64 

A ( A ) = 

0 1 _ 

\ 3 + 

p 591.229 

- 618.051 

• 

a 2 + 

6568.29 

- 5763.7 


[9481.76 


+ 


-8521.031 


8095.49 


7263.08 
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The solvents for A(A) are 


R, 


Rn 


R„ 


1 -2 
3 -4 


det[M-R^] - 


-15 6 

-5 -2 


det [ AX-R 2 ] « 


360 

-480 

320 

-424 


det[M-R 3 ] 


X 2 +3A+2 - (A+l) (A+2) 
A 2 +17X+60 - (A+5) (A+12) 

- A 2 +64A+960 - (A+24) (A+40) 


The system equation x(t) * Ax(t)+Bu(t) can be solved by considering the re' 
duced equations 


(7.14) v^t) - R lV (t) + B^t) 


provided that the overhead for computing the Krylov transformation and that 
for finding the solvents are acceptable as to efficiency and accuracy. 


8. Summary 


The mathematical analysis in this report had two objectives: the first 

to bring together the mathematical tools for understanding matrix poly- 
nomials, with the second of applying these tools to decoupling of system 
equations. The study of matrix polynomials, or lambda matrices, is justified 
as vibrating systems are generally defined with second-order matrix poly- 
nomials, see [19]. Although there are several books and journal articles 
on matrix polynomials, a library search did not reveal a single source of 
the material that is complete. It is essential that the mathematics of 
matrix polynomials be understood before efficient algorithms for analyzing 
vibrating systems can be developed. Present algorithms are not capable of 
handling large space systems — systems with the number of modes greater than 
1000. 

The report is incomplete in several areas. The algorithm for determin- 
ing latent projectors from latent vectors for matrix polynomials with repeated 
roots was not fully developed. Algorithms for efficient computation of 
latent roots and latent vectors have not been devised. A thorough literature 
search has not revealed the availability of a computer program for that 
purpose. In addition, the decoupling scheme was presented but the computa- 
tional algorithm for that task was not described. Work will continue in 
these areas. 

The damping of large-space structures is an important engineering design 
task that must be addressed and algorithms must be developed for that purpose. 
The analysis of matrix polynomials as well as a comprehensive understanding 
of how damping affects the overall mathematical structure is an integral 
part of designing large space structures. Modification of matrix polynomials. 
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which is necessary for the Inclusion of damping, was not considered In the 
report. Studies will begin In that area with the emphasis on second order 
polynomials. 

This report Is only a beginning for the several tasks described above. 
The development of algorithms for analysis, and design of large space 
structures will be addressed during the next year of work. 
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Appendix Krylov Trans formation 


The Krylov transformation Is a useful algorithm that transforms a 
general mnxmn matrix A to the companion form. Let denote the companion 
form, then 


(A.l) 


A * K A K 
c 


-1 


The Krylov transformation Is a similarity transformation that leaves the 
eigenvalues invariant but changes the eigenvectors of A from Q to Q c , where 


(A. 2) 


and 


(A. 3) 


Q “ 


Q 11 A 1 


q C = 


■ Q u 

Q 12 

• • » 

I 

A 

or 

<»21 

• 

q 22 

• 

• • • 

• • • 

q 2m 

• 

• 

_ q ml 

• 

q m2 

• • • 

ii • • 

• 

Q 

mm J 

f Q ll 


q !2 

• 


Q 12 A 2 


0 A 1 ” - ^ 0 

L q il A l 9 12 A 2 


m 


°lm\n 


Q, A 11 *" 1 
lm m 


when A has distinct eigenvalues. The similarity transformation in (A.l) and 
the structure of (A. 2) and (A. 3) requires that K satisfy 


(A. 4) 


Q c - KQ 
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Now from (A. 2) - (A. 4), if m - 3 



‘ou 

Q 12 

Q 13 


1 

H* 

K 12 

5 * 

b* 

U> 

l 



Q 12 

Q 13 

(A.5) 

Q 11 A 1 

Q 12 A 2 

Q 13 A 3 

M 

K 21 

K 22 

K 23 



Q 22 

Q 23 


^ Q 11 A 1 

Q 12 A 2 

Q 13 A 3 _ 


, K 31 

K 32 

K 33. 


Q 31 

m 

Q 32 

Q 33. 


The nxn matrix Q^A.^ is from (A. 5) 


Q 11 A 1 “ K 21 Q 11 +K 22 Q 21 +K 23 Q 31 


or 


(A. 7) 


Q 11 A 1 Q 11 “ K 21 +K 22 Q 21 Q 11 +K 23 Q 31 Q 11 " A B1 


where A^ is the matrix that would be in the upper diagonal from the simi- 
larity transformation T AT given in Section 7. The matrices and 

^31^11 are “ R i2 an£ * ^13 as defined Section 7 but for the general matrix 
with 


(A. 8) 


T 


( 2 )' 


X 

CN 

r-t 

erf 

i 

R 13 


I 

“ Q 12 Q 22 

Q 13 Q 33 

R 21 

-i 

R 23 

1 

~ 4 

^21^11 

-I 

Q 23 Q 33 

_ R 31 

" R 32 

I 


. Q 3i Q n 

-Q32Q22 

I 


Analysis of the other equations in (A. 5) gives the set of equations 


A B1 K 21 +K 22 R 21 +K 23 R 31 
*12*8 2 = K 21 R 12 +K 22 +K 23 R 32 


(A. 9) 


R 13 A B3 ~ K 21 R 13 +K 22 R 23 +K 23 


ti 


(A. 10) 


1 " K 11 +K 12 R 21 +K 13 R 31 
R 12 * K ll R 12 +K 12' fK 13 R 31 
R 13 “ K 11 R 13 +K 22 R 23 +K 23 


A B1 1,5 K 31 +K 32 R 21 +K 33 R 31 
(A. 11) R i2 A B2 “ K 31 R 12 +K 32 +K 33 R 32 
R 13 A B3 “ K 31 R 13 +K 32 R 23 +K 33 


The Krylov matrix is given from (A. 9)- (A. 11) as 



*11 

K 12 

*13 ! 

' (A. 12) 

K 21 

k 22 

k 23 ': 


_ K 31 

K 32 

K 33 j 

or 

K 11 

*12 

K 1 3 1 

(A. 13) 

K 21 

K 22 

1 

K 23 


/31 

K 32 

K 33_ 


R 12 R 13 


21 


R, 


2.3 


R 31 R 31 1 


R. 


12 


B1 

2 

Bl 


R 


13 


A B1 R 12 A B2 R 13 A B3 
L A B1 R 12 A B2 R 13 A B3, 


R 12 

D A 
K 12 rt B2 


R 13 

R 13 A B3 


R 12 A B2 R 13 A B; 


R 12 r i; 


21 


L R 31 R 32 


R 2 

1 J 


'll 


12 


13 


L A 11 +A 12 A 21 +A 13 A 31 A 11 A 12 +A 12 A 22 +A 13 A 32 A ll A 13 +A 12 A 23 +A 13 A 33 


V 1 

e r A 

k A 2 


= K 


I 

] 

J 
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where * [10 0]. The similarity transformation given in (A.l) is required 
to show that the Krylov matrix is as given in (A. 13). 

The general form of (A. 13) is 


(A. 14) 

I 

l 

I 

I 

I 

| where A is the general matrix to be reduced to the companion form, A . 

l 

1 The inverse of all matrices must exist. Row-column interchanges of 

A can be made in most cases to assure the existence of A . The algorithm 

C 

will not have good accuracy wh< n m is large particularly when the matrix 
is stiff} the eigenvalues have a large spread in magnitude. 
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